The nitrogen-vacancy center in diamond re-visited 
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Symmetry considerations are used in presenting a model of the electronic structure and the asso- 
ciated dynamics of the nitrogen- vacancy center in diamond. The model accounts for the occurrence 
of optically induced spin polarization, for the change of emission level with spin polarization and for 
new measurements of transient emission. The rate constants given are in variance to those reported 
previously. 
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I. INTRODUCTION 

Since the nitrogen-vacancy (N-V) center in diamond 
has been detected at a single site level 0, 0, the cen- 
ter has attracted attention for various quantum informa- 
tion processing applications. For example, the center has 
been used as a single photon source for quantum cryp- 
tography 0, Q and work in this area has included im- 
pressive demonstrations Q. Another area of interest re- 
sults from the center having a non-zero spin ground state. 
The ground state spin can be the qubit and the optical 
transitions utilized for readout and for qubit manipula- 
tion in quantum computing applications. There are a gain 
impressive demonstrations in this area 0, IS IS EH El • 
With these successes and additional programs under de- 
velopment it might be expected that the properties of the 
N-V center are well understood. However, this is not the 
case. Despite the extensive publications the center is not 
well understood and the literature contains many inaccu- 
racies. The purpose of this paper is to re- visit our knowl- 
edge of the nitrogen-vacancy center, provide an account 
of the electronic energy levels and explain the dynamics 
of the center under optical excitation. 



II. NITROGEN- VACANCY CENTER 

The N-V center occurs in diamond containing sin- 
gle substitutional nitrogen when irradiated and annealed 
|l2l H^l . Electron irradiation with energies greater than 
200 keV creates vacancies EH ■ The vacancies are mobile 
at 800° C and can become trapped adjacent to the nitro- 
gen impurities. The nitrogen- vacancy complex formed 
has a strong optical transition with a zero-phonon line 
at 1.945 eV (637nm) accompanied by a vibronic band at 
higher energy in absorption and lower energy in emis- 
sion. The zero-phonon line has been studied by Davies 
and Hamer E3- They have studied the effect of uni- 
axial stress and from the splitting and polarization they 
established that the transition is associated with an or- 
bital A - E transition at a site of trigonal symmetry. 
The trigonal symmetry is consistent with an adjacent 
nitrogen- vacancy pair with C^ v symmetry. In other stud- 
ies using optical excitation Loubser and van Wyk de- 



tected electron spin resonance (ESR) signals of a spin 
polarized triplet (S = 1) p|. The ESR they observed 
was associated with a center having trigonal symmetry 
and the magnitude of the optically induced signal was 
found to vary with wavelength in correspondence with 
the strength of the A - E optical transition. Hence, the 
center was attributed to the same nitrogen- vacancy com- 
plex. With an integer spin (S = 1) the center must have 
an even number of electrons and it is taken that the neu- 
tral nitrogen- vacancy complex with five electrons has ac- 
quired an additional electron from elsewhere in the lat- 
tice, probably from another substitutional nitrogen atom. 
There will then be six electrons occupying the dangling 
bonds of the vacancy complex 15]. Loubser and van 
Wyk proposed that the spin polarization arises from a 
singlet electronic system with inter-system crossing to a 
spin level of a meta-stable triplet. However, it was estab- 
lished from hole burni ng Jlol , optically detected magnetic 
resonance fiffl, ESR [l£| and Raman heterodyne mea- 
surements 19] that the triplet is the ground state rather 
than a meta-stable state. Therefore, their model has to 
be modified to give a 3 A ground state and a 3 A ~ 3 E op- 
tical transition. The optically induced spin polarization 
can still arise from inter-system crossing and an account 
is given in this paper. 

The six electrons occupy the dangling bonds associ- 
ated with the vacancy complex. A discussion of this sit- 
uation is included in a treatment by Lenef et.al. l20l . 
Although they did not adopt the six electron model [21| 
they did give a very useful general treatment including 
the six electron situation and their presentation allows 
the present discussion to be brief and more descriptive. 
The dangling bonds are formed from sp 3 orbitals of the 
carbon and nitrogen atoms and in a vacancy approximat- 
ing Td symmetry these can be combined to form a\ and 
t2 molecular orbitals with A± and T2 symmetry, respec- 
tively 22]. From symmetry and charge overlap consid- 
erations the a\ is considered to be lower in energy and 
the £2 higher. With six electrons a\t\ will be the lowest 
energy configuration and this can also be described as 
a t\ hole system. When one of the neighboring carbons 
is replaced by a nitrogen the T^ symmetry will be low- 
ered to trigonal and the ti orbital will be split to give, in 
C-iv notation, an a\ and e orbital. The e hole is lower in 
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energy and, hence, the lowest energy configuration will 
be e 2 , next lowest ea\ and the a\ highest. The spin- 
orbit wave functions for the e 2 configuration give 3 ^2, 
1 Ai and 1 E states and the ea\ configuration 3 E and 1 E 
states. The a\ gives an 1 Ai. The optical transition is 
associated with triplets and so the ground state is at- 
tributed to the 3 A2(e 2 ) state and the excited state to the 
3 E(eai) state. There are singlets 1 Ai(e 2 ) and 1 £'(e 2 ) and 
1 E(aie) which could lie in the same energy range as the 
triplets. The 1 A 1 (af) lies higher. It has been assumed 
that the 1 A 1 (e 2 ) lies between the triplets and the present 
treatment accepts this assertion and will be shown to be 
consistent with observation. The possibility of interme- 
diate 1 E level(s) will be discussed. 

The energy, V, of each state 3 ^2, 3 E, 1 A\ is determined 
by the above bonding considerations and the Hamilto- 
nian including spin-orbit, Vso, an d spin-spin, V55, in- 
teraction is given by: 

H = V + V so + V ss (1) 

Spin-orbit and spin-spin do not affect the degeneracy 
of the singlets whereas the spin degenerate ground state 
is only affected by spin-spin interaction normally written 
as: 

Vss = pS 2 z +p(Sl + S 2 y ) (2) 

where p and p are the axial and non-axial coeffi- 
cients. The interaction spits the ground state into a sin- 
glet, |A2,S Z > with symmetry Ai, and doublet, |A2,S K >, 
A2,Sj,> with symmetry E. The spin states |S 2 > and 
(|S a >, |S V >) are not mixed. 

The 3 E state are affected by both terms. Vso has the 
form: 

V so = X(L Z S Z ) + \'(L X S X + L y S y ) (3) 

where A and A are the coefficients associated with the 
axial and non-axial spin-orbit interactions. The axial 
A(L Z S Z ) spin-orbit interaction splits the 3 E spin triplet 
into three two-fold degenerate states; E, E', and an (Ai, 
A2) pair. Within the 3 E state non-axial spin-orbit in- 
teraction is small and will be neglected at present. As 
will be discussed shortly spin-orbit can give mixing be- 
tween adjacent states. This can cause a shifting of levels 
which can be calculated by including this state in the cal- 
culation or by including spin-orbit interaction to second 
order. The form of this latter term is the same as that 
of spin-spin interaction. In high symmetry (eg T^) this 
has the form p^ : 

Vss = p[{L.Sf + 1/2{L.S) + L(L + l)S(S + I)] (4) 

In trigonal symmetry one has to allow for the differ- 
ence in the axial and non- axial terms. The interaction 



modifies the separation of the above states but does not 
change the wave- functions and the wave functions are the 
main interest here. The states and the symmetry adapted 
wave functions are given in Fig. I (a). The states with S z 
spin projection are not mixed with the spin states with 
S x and S y spin projection. The 3 A2 «-> 3 E transition is 
orbitally allowed and, as spin projection is not changed 
by the electric dipole operator, the optical transitions 
will be the same strength for each of the spin projections 
(shown as solid vertical arrows in Fig. 1(a)). There are no 
transitions involving a change of spin and it can be con- 
cluded that optical cycling of the 3 A2 - 3 E transition will 
not result in change of spin projection and consequently 
can not give any spin polarization. This situation does 
not change when considering vibronic interactions as spin 
projection is conserved. 

Spin-orbit interaction mixes singlets and triplets which 
transform according to the same irreducible representa- 
tion. The mixing provides an avenue whereby symmet- 
ric vibration can cause a population relaxation between 
the mixed singlets and triplets. The symmetry consid- 
erations, therefore, determine the allowed inter-system 
crossing and these are also shown in Fig. 1(a). Where 
there is only a x Ai singlet level the inter-system crossing 
is restricted to states that transform as Ai irreducible 
representations. There can be excitation of population 
out of the E ground state level with S x or S y spin projec- 
tion to the Ai spin-orbit level of the 3 E state. Population 
in this state can decay via the singlet to the S z spin pro- 
jection of the ground state. Such an excitation and de- 
cay, therefore, causes a re-orientation of the ground state 
spin projection and with continuous excitation popula- 
tion can be transferred to the ground S z spin state. This 
is consistent with the preferred spin orientation estab- 
lished experimentally [2J|. With these selection rules, 
assuming the optical induced spin polarization is faster 
than spin-lattice relaxation, the spin polarization would 
be 100% but this is not what is observed |25|. 

The extra consideration is the non- axial A (L^S^ + 
LySj,) spin-orbit interaction which we have previously 
neglected. This spin-orbit term causes a mixing of the 
two 3 E states (denoted E and E') transforming as an 
E irreducible representation. These states have differ- 
ent S z and S x ,S y spin projections and, hence, the mixing 
gives rise to optical transitions that do not conserve spin 
(dashed arrows in Fig.l}. The transitions are observed in 
hole burning spectra Ha |H [H III El E EI and are 
the transitions that limit the degree of spin polarization. 
It is anticipated that the strength of the non-axial spin- 
orbit interaction is small. This is concluded from con- 
sideration of the the one electron operators. Spin-orbit 
interaction will be isotropic for t2 orbits in T^ symmetry 
but the axial contribution is quenched when the orbit is 
split into two- fold degenerate and non-degenerate states. 
In the non-axial spin-orbit interaction between the 
individual two-fold degenerate states is formally allowed 
but its contribution arises from higher order effects. The 
non-axial spin-orbit interaction is, therefore, anticipated 
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to be small and the non-spin conserving optical transi- 
tions will be weak. Although weak the transitions play an 
important role in limiting the degree of spin polarization. 



A1 A1 a |E X , S„>*|Ey, Sy>\ 
A2 A2 = |E X . S y >- Ey, Sx>/ 




E Ey =|A 2 ,S X > 
Ex=|A2 ,S y > 

A1 Ai =|A2 , S 2 s 



FIG. 1: (a) Energy levels of the N-V center in C3,, sym- 
metry. The excited state is split by spin-orbit interaction 
whereas the ground state is split by spin-spin interaction (not 
shown to scale). The figure gives the symmetry adapted wave- 
functions. The solid arrows indicate the spin-allowed optical 
transitions. The dashed arrows indicated weak transitions 
which are allowed through the mixing of the ( 3 E)Eand ( 3 £)E' 
basis states by transverse spin-orbit interaction. The diago- 
nal arrows give inter-system crossing allowed by spin-orbit 
interaction, (b) Energy levels and wave functions of the N-V 
center in the presence of a strain held. The wave functions 
are appropriate for a strain field perpendicular to a reflection 
plane. The transitions are derived from those allowed in (a). 

We briefly consider the situation where a 1 E state(s) 
lies between the 3 A2 and 3 E. If this were the case there 
could be symmetry allowed relaxation from the excited 
triplet ( 3 E)E and ( 3 E)E' states to such an intermedi- 
ate *E singlet state. However, using the two hole de- 
scription the relaxation is found to be forbidden in the 
case of the triplet ( 3 E)E' state (cancellation of terms for 
a 1 E(aie) and forbidden for one electron operators for 
1 E(e 2 ) states) implying that there will still not be pop- 
ulation transfer out of states with S z spin projection. 
There can be relaxation from the other triplet ( 3 E)E 
state involving the S x , Sj, spin projections to a 1 E level. 
Given that there is also an intermediate : Ai there can 
be two alternate situations depending on the ordering of 
the 1 Ai and 1 E singlet levels. Should the : E state lie 
lowest the decay will be to the ( 3 A2)E component of the 
ground state. This state has a Sj, S y spin projection 
and the optical cycling will have had no consequence im- 
plying no change in spin orientation. However, if the 
1 E state lies above the 1 A\ level there will be radiative 
(or non-radiative) decay to the lower singlet level, the 
1 A\ 1 followed by the relaxation, as discussed in the pre- 
vious paragraph, to the S z spin projection of the ground 
state. This process will lead to the same spin change and 
spin polarization as before. The simple consequence of 



involving the : E singlet is that the total ( 1 ^4i plus 1 E) 
spin polarization process will be more efficient. Thus, if 
there is an intermediate 1 E state, to be consistent with 
observation it must lie at an energy higher than the 1 A\ 
state. It is recognized that this order is in disagreement 
with calculated energy levels |22j. However, there is an 
appeal of including a higher single 1 E state as it could 
lie close to the excited triplet level and the 1 Ai close to 
the ground state thus accounting for the relatively fast 
inter-system crossing reported below. However, there is 
no fundamental difference in the dynamics and it is suffi- 
cient in this work to restrict the discussion to one singlet 
level, a : Ai. 

It is common for there to be strain in diamond and 
it is worth considering the consequence to the energy 
levels and the associated dynamics. There will be no 
fundamental change with axial strain as all it gives is a 
uniform shift of the energy levels and no change of the 
wave functions or the dynamics. However, the compo- 
nent of the strain at right angles to the axis of the center 
lowers the symmetry and the extra crystal field lifts the 
orbital degeneracy of the excited state to give two orbital 
non-degenerate states denoted E x and E y [2(J. We con- 
sider the case where this strain splitting is larger than 
the spin-orbit interaction. Where the strain retains a re- 
flection plane the Z, X, Y axes will be determined by 
symmetry. The wave functions for this situation are as 
given in Fig. 1(b). The diagonal spin-orbit interaction is 
quenched and the order of states (same in both optical 
components) are determined by spin-spin interaction. In 
the limit of small non-axial spin-orbit interaction there 
is still little mixing of the S^, S y spin states with the 
S z spin states. The transitions and inter-system crossing 
can be determined from the parent state and are shown 
in Fig. 1(b). The spin allowed transitions will have near 
the same strength as in the zero strain case but they are 
now totally polarized. The only minor change is loss to 
the dashed transitions between the states with different 
S x and S y projections. These arise where the strain is 
not sufficient to totally quench the effect of the diagonal 
spin-orbit interaction. A more significant effect is in the 
strength of the non-spin-conserving transitions induced 
by non-axial spin-orbit interaction. They will become 
significantly stronger as the separation of the S z and 
(S x , S y ) is reduced and the mixing increased. When the 
strain does not retain reflection symmetry the effective X 
and Y axis may be different between ground and excited 
states and all transitions and inter-system crossings will 
become allowed in principle (no symmetry restrictions). 
However, the selection rules will be dominated by those 
allowed in zero order and the perturbation approach in 
Fig. 1(b) is anticipated to give a reasonable approxima- 
tion to the dominant excitation and decay channels. 

Which of the diagrams, Fig. 1(a) or Fig. 1(b) (or an in- 
termediate case), is appropriate for a given center de- 
pends on the relative magnitude of the stress, spin-orbit 
and spin-spin interaction. Spin-orbit for the carbon atom 
is known to be a few cm -1 (~200 GHz) and a spin-orbit 
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splitting of 1 cm -1 (30 GHz) was obtained for the 3 E 
state from optical magnetic circular dichroism measure- 
ments |l (tI | . This value maybe marginally high as an op- 
tical line width of 15 GHz has been reported recently 
for small ensembles within single crystals of diamond 
pflj and there is an indication that spin-spin interaction 
may also contribute to the line width. Santouri et. al. 
|29j have also shown, using two-laser hole burning exper- 
iments, that for even a small strain splitting of 10 GHz 
the energy scheme is equivalent to Fig. 1(b). As strain 
splitting is usually much larger (> 100GHz) it is taken 
that this energy scheme will be typical of centers in an 
ensemble sample as used here. 

The dynamics associated with optical excitation can 
be determined without detailed knowledge of the excited 
state energy levels. This is because from the above dis- 
cussion it can be seen that the system can be reasonably 
quantized by its spin projection S z or (S^, S y ) and this 
is not altered by stress. The dynamics are largely de- 
termined by the crossing between these two spin projec- 
tions and the symmetry considerations give the two main 
mechanisms. The changes of spin are through the inter- 
system crossing via the singlet and through the weak non- 
spin-conserving optical transitions. The effective energy 
scheme is given in Fig. 2 where the S z states are shown 
on the left and the (S x , S y ) on the right. The singlets are 
drawn centrally. Spin polarization involves displacement 
of population from the states on the right to the states 
on the left. The transitions including those that change 
the spin projection are shown as solid lines as determined 
by the parent states in C3„ symmetry. For completeness 
other transitions allowed in low symmetry are shown as 
dashed lines. The low symmetry situation is similar to 
that proposed by others pAj . l33l.l34j and to assist compar- 
ison we adopt their shortened notation. The ground spin 
states arc defined as x, y and z; the excited triplet states 
as x', y' and z' and the singlet level as s. In the following 
experiments we determine the values of the parameters 
using bulk samples. 



III. EXPERIMENTAL DETAILS 

Several type lb diamonds were used in these studies. 
They were irradiated with energetic electrons (> MeV) 
and annealed. The nitrogen-vacancy concentrations var- 
ied from 3 x 10 18 per cm 3 to 10 17 per cm 3 and there 
was no indication that the dynamics of the optical cycle 
varied significantly for concentrations in this range. 

The experiments involve transient and CW excitation 
studies using a frequency doubled Nd:YAG laser. The 
excitation wavelength is 532 nm which is near the peak 
of the 3 A2 <-> 3 E absorption band. The emission was 
detected by a S-20 photomultiplier or pin diode with a 
response time of 30 ns. With the exception of the mea- 
surement of the spectrum the emission wavelength de- 
tected was selected using absorptive filters. 

The studies involved three different excitation and de- 




FIG. 2: Energy levels for a perturbed N-V center. The sig- 
nificant energy levels of Fig. 1(b) are re-drawn to highlight 
the inter-system crossing. All the radiative and non-radiative 
transitions are shown. The allowed transitions considered in 
the text are shown by solid arrows. The dashed vertical ar- 
rows on the right are allowed but can be taken to be zero for 
present experiments. Also the inter-system crossing indicated 
by the dashed arrows are zero in first order. 

tection geometries. For low level excitation the intensity 
was determined from the laser power and the beam diam- 
eter. The emission at 45° was focussed on the pin diode. 
For high intensity measurements a confocal arrangement 
with an oil immersion objective was used. The excited 
diameter was 0.2 micron with the emission spot focussed 
onto a pinhole. A 1 mW beam gave an estimated power 
density of 2 x 10 6 W/cm 2 . With the small excitation and 
detection volumes the signals were weak and long collec- 
tion times were required to obtain satisfactory signal to 
noise ratios. The majority of the measurements were 
taken with a third geometry where satisfactory signals 
could be obtained more rapidly. The light was focused 
with a microscope objective and the back emission col- 
lected by the same objective. For the same laser power 
the intensity at the sample was approximately two orders 
of magnitude lower than that obtained with oil immer- 
sion objective and it required a 100 mW laser beam to 
obtain a maximum intensity of 2 x 10 6 W /cm 2 . The un- 
satisfactory feature of this geometry is that the intensity 
was not constant over the collected volume. 

At low intensities where slow (>ms) responses were 
obtained the light was gated with a mechanical chopper 
whereas for the fast speeds associated with the high in- 
tensities the light was gated using two acousto-optic mod- 
ulators in series. The rise time of the A-0 modulators is 
30 ns. When gating off there was a weak component that 
lasted for a /is. However, no measurements were taken 
when gating the light off. 
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A. Preliminary measurements 



The emission of the N-V center has been reported on 
many occasions. The general characteristics are shown 
in Fig. 3. Near room temperature the emission gives 
a band stretching from 630 nm to 800 nm with vibra- 
tional structure and a weak zero-phonon line at 637 nm 
(Fig. 3 (a)). Cooling has little effect on the vibronic ab- 
sorption band |3^| and consequently there is little change 
in the amount of light absorbed from the laser beam when 
using an excitation wavelength within the vibronic band. 
Consistent with this, the total emission does not show a 
significant change when the temperature is lowered (in- 
sert in Fig. 3(a)). The most obvious change in cooling to 
low temperatures is that the zero-phonon line becomes 
sharper and more prominent (Fig. 3(b)). The Huang- 
Rys factor is 3.7 with only 2.7 percent of the transition 
strength being associated with the zero-phonon transi- 
tion [35| . At high excitation densities there can be photo- 
ionization of the N-V center of interest and the creation 
of the neutrally charged [N-V]° center [3(|[37|]. This cen- 
ter has a zero-phonon line at 575 nm with a vibronic band 
to lower energy 38] . The increased contribution of this 
center at high intensities is illustrated in Fig.3(b). 

As discussed earlier the N-V center exhibits optically 
induced spin polarization of the ground state triplet and 
there is a change of the emission intensity associated with 
this polarization (no change in absorption). The change 
is illustrated in Fig. 4. In the dark the sample becomes 
unpolarized and when excited the emission has an initial 
intensity level A. After exciting for a period the sample 
becomes polarized and the strength of the emission in- 
creases to a second value B (Fig. 4(a)). The rate at which 
the emission increases from the level A to B is linearly 
dependent on the excitation intensity. Should the sample 
be in the dark for a period less that the spin-lattice relax- 
ation time (Ti) the initial emission level will be at a value 
between A and B and varying the dark period varies this 
level. This variation of emission level with the duration 
in the dark can be used to establish the spin lattice relax- 
ation time and measurements of this type are shown in 
Fig. 4(b). The A/B ratio is of considerable interest. To 
ensure a satisfactory A measurement it is necessary to 
have the sample in the dark for a period long compared 
to Ti whereas to obtain the saturation value B the in- 
tensity has to be sufficient to obtain the higher emission 
level within a time short compared to Ti. A decrease of 
the spin polarization through spin diffusion must also be 
avoided. When meeting these conditions the A/B ratio 
is 0.86 ±0.02. 



B. Dual pulse measurements 

The sample is excited with two intense excitation 
pulses and the delay between two pulses is varied. These 
measurements utilized the confocal arrangement with oil 
immersion and intensities of 10 6 W/cm 2 . The repetition 
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FIG. 3: Spectral characteristics of the N-V Center, (a) Emis- 
sion spectra of a diamond containing a high concentration of 
N-V centers measured near room temperature using low ex- 
citation densities of 1 W /cm 3 . Insert: the variation of the 
N-V emission at 666 nm as a function of sample temperature 
(b) Low temperature emission spectra of the N-V center mea- 
sured with various excitation densities (from reference |3(|). 
Insert: variation of N-V emission intensity as a function of ex- 
citation energy density (units of 1 x 10 5 W /cm 3 ). Emission 
normalized to excitation intensity. 



rate was 10 kHz. The results are shown in Fig. 5(a). In 
the figure the response associated with the first pulse of 
every pair is overlapped whereas the emission associated 
with the second pulse is displaced as the delay between 
the pulses is varied. Every pulse exhibits an initial peak 
followed by a decay within a /js to a lower level. For 
the first pulse the magnitude of the peak is 1 /3 that of 
the equilibrium signal. When the light is gated off the 
sample has to remain in the dark for a period before the 
second pulse gives a peak and the magnitude of the peak 
increases as the duration of the dark period is increased. 
This recovery has two components and the faster recov- 
ery is shown to have a response time of 0.3 /is (Fig. 5(b)). 
There is a slower recovery over the 100/is between pulses 
and we will comment on this slower recovery later. 

For the second series of measurements the emission is 
obtained for a pulse pair of 1/zs duration separated by 
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FIG. 4: Emission obtained when gating on excitation at time 
0. The excitation is at 532 nm and is gated with a mechanical 
shutter, (a) Response for various laser powers; 10 mW corre- 
sponds to an intensity of 1 W/cm . The emission increases to 
an initial value A, 86% of its final value. The rate of increase 
of the emission to its final CW value B depends on intensity 
of excitation. Sample has spin-lattice relaxation time of 100 
ms and sample held in dark for 500 ms prior to excitation, (b) 
Responses for various periods in the dark. For these traces 
the sample is cooled and the spin-lattice relaxation time is 
increased to 500 ms. 



but with long delays between the pulse pairs. The 
measurements were made using the alternate geometry 
where the excitation and detection involved a larger but 
less well defined volume. The period between pairs (>10 
ms) was chosen to be much larger than ground state spin- 
lattice relaxation time, Ti (1 ms). The intention is for 
the system to be unpolarized at the start of the first pulse 
but polarized at the start of the second pulse. Fig. 6(a) 
shows the results of the emission response for pulse pairs 
for laser powers from 5 mW to 200 mW corresponding to 
estimated intensities of 10 5 W/cm 2 to 4 x 10 6 W/cm 2 . 
The emission is restricted to longer wavelengths (> 700 
nm) to avoid including emission from [N-V] centers. The 
emission of the [N-V]° center by itself was also obtained 
by detecting the emission at 590 nm using a narrow band 
filter. This emission is shown in Fig. 6(b). 

The responses in Fig. 6 (a) show peaks at the start of 



FIG. 5: Emission response to a double light pulse (a) The 
emission is shown with various delays between the two pulses. 
The responses for the first pulses all overlap whereas the sec- 
ond pulse is delayed by the dark period. The peak associated 
with the second pulse recovers with the dark period and the 
height of this pulse is plotted in (b) as a function of the du- 
ration in the dark. The dashed curve illustrates the response 
for a 0.3 /its recovery rate. The repetition rate of the pulse 
pair was 10 kHz. The intensity of the excitation was 3 x 10 6 
W/cm 2 . 



each pulse and the magnitude increases with intensity. 
The peak height of the second pulse is consistently several 
percent lower than that of the first. Also the decay rates 
arc different for the two pulses, the first being faster than 
the second. 

The effect of applying a weak magnetic field of a few 
hundred gauss was also recorded. The response to a pair 
of intense (100 mW) excitation pulses at 532 nm were 
measured both with and without the magnetic field ap- 
plied in a random direction (Fig. 7). In comparing the 
two traces there are three significant differences. The 
first emission peak is almost the same for the two traces 
but the subsequent decay is to a much lower level when 
the magnetic field is applied. With the second pulse there 
is a significant difference in the peak heights, being lower 
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FIG. 6: Emission of the N-V center for a pair of square wave 
excitation pulse of light at 532 nm. The light is focused with 
a microscope objective and for 100 mW the intensity is 2 x 
10 6 W/cm 2 The back scattered emission is detected using a 
(a) 700 nm long pass filter, (b) 100 nm band pass filter at 590 
nm. 



FIG. 7: Emission of the N-V center for a pair of square wave 
excitation pulses of light at 532 nm with laser intensity of 100 
mW. No magnetic field is applied in the case of the upper 
trace and corresponds to the 100 mW trace in Fig. 6. For 
the lower trace a field of 500 gauss is applied in a random 
direction (not aligned with an axis of any center). 



with the field is applied. Also with the second pulse the 
decay to the lower level is faster when the magnetic field 
is applied. 



IV. RELAXATION AND INTER-SYSTEM 
DECAY RATES 

The dynamics scale with the lifetime of the excited 
state and various values have been reported in the lit- 
erature. Collins et al [3!j have obtained values of 12.9 
± 0.1 ns for a natural diamond and 11.6 ± 0.1 ns for a 
synthetic diamond. Lenef et.al. [20j has also obtained a 
measure of 12.96 ± 0.14 ns obtained in relation to pho- 
ton echo measurements. A value of r = 13 ns is adopted 
here. The radiative rate constants are, hence, k x >x = 
= k z ' z — 77 x 10 6 s _1 . (The present experiments are not 
sensitive to the dashed vertical transitions in Fig. 2 be- 
tween the x, x' and the y, y' states and the associated 
rate constants can be taken to be zero, ie k xy < = k yx i = 
0) . It is noted that the present model indicates that there 



should be two components to the emission decay associ- 
ated with excited states, one with z' and one with x',y'. 
The emission from z' excited state has the slower decay 
rate and dominates when the system is spin polarized, 
perhaps explaining why the two values (differing by only 
30% ) have not been detected. Two components of 9 ns 
and 2 ns have been observed by Hanzana et al |4jj but 
the measurements are inconsistent with previous values. 
These lifetimes require further investigation. 

The increase in emission from a level A to a level B 
can be used to establish the fraction of the population 
transferring into the singlet. For example the 14% emis- 
sion change between having all the population in the z 
spin ground state and having the population evenly dis- 
tributed between the three spin states implies 27 % of 
the population from each of the x', y' excited spin states 
transferring non-radiatively to the singlet, s. This has 
to be increased to 39% to allow for the incomplete spin 
polarization reported below. The inter-system crossing 
rates are then k x / s = k. y < s = 30 x 10 6 s _1 . In correspon- 
dence with the model, the inter-system crossing from the 



8 



z' state to the singlet is taken to be zero and, hence, k z > s 
= 0. 

For the low intensities no population is maintained 
in the singlet. This is changed at high intensities and 
the transient emission displays different characteristics. 
With population being stored in the singlet level there 
is a drop in emission and this is observed in all of the 
two-pulse experiments (Figs. 5, 6 and 7). There will be 
no initial peak if the excitation is gated on and off within 
a few ns as there will be no change in the singlet popu- 
lation. The recovery of the peak requires a period in the 
dark (Fig. 5) and the time required corresponds to the 
rate at which population decays from the singlet to the 
ground state. The value of the singlet lifetime obtained 
from this peak recovery is 0.3 /is and this gives k sz = 3.3 
x 10 6 s _1 . As in the model we take k sx = k sy = 0. 

In a recent paper we have reported that the maximum 
spin polarization obtained for an ensemble under contin- 
uous excitation is 80% |25j. This means that the prob- 
ability of optically transferring spin projection from the 
z spin state to an x, y spin is 1/4 of the above process 
giving rise to the spin polarization. The mechanism is 
attributed to the non-spin conserving optical transitions 
(diagonal arrows between triplet levels in Fig. 2) and 
implies that the rate constants are k zx i = k zy i = k x i z 
= kj,/ 2 = 1.5 x 10 6 s _1 . Loss of spin polarization could 
also arise from the reverse inter-system crossing via the 
singlet level. However, the model predicts that the rate 
is zero, k z i s = and, as given previously, the decay from 
the singlet to the x and y ground states are also zero. 



V. RATE EQUATIONS 

In the previous section estimates of the parameters of 
the model in Fig. 2 have been obtained from simple ex- 
perimental observations. By adopting these parameters 
we can determine the populations and the emission for 
any optical field by solving the classical rate equations: 



drill dt = ^j(kji n j ~ kijUi) (5) 

where n^ is the population of level i (i = z, x, y, z', 
x', y', s). and k^ gives the rate for the i — ► j transition. 
The significant parameters associated with the optical 
transitions are: 



[kz'zi k x ' X: fey'y)i \]^x'z% ky'z: k z > Xl kz'zyji (k x 'y : ky'x) (6) 

where the values within brackets are equal in first or- 
der. The related optically driven terms are obtained by 
setting kij — k* kji where k indicates the strength of the 
optical pumping, k = 1 corresponds to the case where 
the optical pumping rate of the allowed transitions equals 
the emission decay rates. The inter-system crossings are 
determined by the rates: 




10000 12X0 14000 



Time |is 

FIG. 8: Emission calculated from solution of the rate equation 
for energy scheme in Fig. 2. The intensities k are given in 
units of 1/r. The value of the parameters in units of 10 6 s _1 
introduced in the text are k zz i — k xx > — k yy , — 77; k sz — 3.3, 



ksx - 
= k 7 . 



0; k z 



0; k. 



yy' 
30; k z 



= 1.5 



(A.v.s), (kx'si ky' s ), (k sz ), (ksx-, k S y) 



(7) 



There are no reverse terms associated with the inter- 
system crossing and, hence, the related parameters with 
the indices reversed are zero. Likewise all the relaxation 
rates between the spin levels z, x and y and between 
levels z', x' and y' are small. These parameters can be 
considered equal and given a small value but the effects 
are not significant in the calculated responses. 

Where a sample has been in the dark for a period long 
compared to the spin-lattice relaxation time the popu- 
lation will initially be equally distributed over the three 
ground state spin levels, z, x and y. Emission is estab- 
lished in a time of the order of the excited state lifetime 
of 10 - 20 ns and with continuing excitation the emission 
level increases as population is transferred to the z state 
(Fig. 8). This behavior is in correspondence with obser- 
vation (Fig. 4). However, little significance can be drawn 
as there has not been an independent measurement of 
the optical pumping rate and the magnitude of the rise 
in Fig. 4 has been used in determining the parameters of 
the system. 

Other than the optical pumping rate there are no 
free parameters when calculating the emission associated 
with the dual pulses experiments. In these experiments 
the intensities are high and the transfer of population 
into the singlet level is initially faster than the relaxation 
from the singlet to the ground state. Consequently there 
is a build up of singlet population and associated with 
this there is a drop in the emission level. The situation 
is calculated for a light field switched on and held con- 
stant for 1/is, switched off for \\x and then switched on 
again for a further 1/is. With the intense excitation the 
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system reaches equilibrium during the 1/xs pulse and so 
is spin polarized well before the end of the first pulse. 
In the dark period there is a relaxation of the singlet 
population but there is no loss in spin-polarization. The 
result is that for the second light pulse the system starts 
spin polarized with a preferential population in the z spin 
ground state. In this case the transfer to the singlet is 
less efficient and the build up of population in the sin- 
glet level is slower. This accounts for the observed slower 
drop in emission intensity with the second pulse. The 
behavior for representative intensities is shown in Fig. 9. 

The results of the calculations can be compared with 
the experimental measurements of Fig. 6. It should be 
recognized that the calculations are for a simpler sit- 
uation than realized experimentally. The calculations 
are for an ensemble of identical centers with identical 
optical pumping rates whereas the experiment involves 
four orientations and variation in the optical pumping 
rates. The consequence of these factors can be approxi- 
mated by adding a square wave emission response to the 
calculated emission response before comparing with ex- 
periment. The structured component of the response ( 
the 'peak') will then represent a smaller fraction of each 
pulse. Another important consideration is that photo- 
ionization has not been included. In the calculations the 
peak of the second pulse is stronger than that of the first 
whereas it is the reverse in the experiment. The differ- 
ence is due to photo-ionization. Photo-ionization causes 
there to be a reduction in the number of N-V centers dur- 
ing the first pulse (giving small alteration to the slope). 
The recovery is slow and there is no recovery during the 
short dark period. Hence, the number of centers involved 
is larger at the start of the first pulse than at the start 
of the second pulse. When allowing for these factors and 
recognizing that the parameters have been determined 
from independent measurements the correspondence be- 
tween the calculated responses in Fig. 9 and equivalent 
experimental traces in Fig. 6 is very satisfactory. 



VI. MAGNETIC FIELD CALCULATION 



A magnetic field other than a long the trigonal axis 
causes mixing of the spin states | l5j and consequently 
with a randomly oriented field the populations are not 
associated with separate z, x and y states. The effect 
can be approximated by retaining equal populations in 
the three spin projections and the result of doing this is 
shown in Fig. 10. The upper trace gives the response 
in the absence of a magnetic field and is the same as in 
the previous section with k = 1. For the lower trace the 
populations in the three ground spin states are equalized. 
When this is done to approximate the effect of a magnetic 
field, there is no spin polarization and the responses are 
the same for the two pulses. With the field applied there 
continues to be optical excitation from the x, y state and 
more efficient transfer into the singlet level. The result 
is that the equilibrium population in the singlet level is 
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FIG. 9: (a) Emission predicted from the solution of the rate 
equations for a pair of excitation pulses. The emission is 
shown for various excitation intensities k shown on the right 
in units of 1/r. The value of the parameters are given in the 
text and summarized in the caption of Fig. 8. 



higher and less centers contribute to the emission. The 
final emission is lower and this is what is observed. It 
should be noted that equivalent effects can be obtained 
by applying resonant microwave fields to maintain pop- 
ulation in the x, y states. At high intensities the drop in 
equilibrium emission level caused by the microwave field 
can be much larger |4l| than the 14% obtained at low 
intensities. This is due to the change of the population 
stored in the singlet level. 

In comparing the responses with and without an ap- 
plied magnetic field there is a variation in the magnitude 
of the peaks. The peak in the emission associated with 
the first pulse is similar with and without field. How- 
ever, the magnitudes associated with the second pulse 
are very different. As discussed previously, the difference 
in peak heights between the first and second pulses is 
due to photo-ionization varying the number of centers. 
A smaller second peak indicates that the magnetic field 
has caused an increase in the photo-ionization. This can 
be attributed to the field increasing the population in 
the excited states (excited triplet plus singlet) and the 
ionization being out of these states. It is desirable to es- 
tablish whether the ionization occurs through tunnelling 
out of these states or is light induced. This requires fur- 
ther investigation. 
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FIG. 10: Emission predicted from solutions of rate equations 
illustrating effect of an applied magnetic field. The solid line 
in the upper trace shows the emission of system determined 
from solution of rate equations for energy scheme as shown in 
Fig. 2. The dashed line indicates the variation in population 
of singlet level. For the lower trace the three ground states are 
mixed, effectively maintaining equal population in the three 
spin projections. The excitation rate is k = 1 in units of 
1/t. The value of the parameters are given in the text and 
summarized in the caption of Fig. 8. 



VII. COMPARISON WITH OTHER WORK 

There have been many publications referring to the 
singlet level in the N-V center and many of these publica- 
tions give information that is in conflict with the current 
model. The disparities are discussed below. 

In the present work it is shown that the singlet has a 
short lifetime of 0.3 fis. This contrasts with the value of 
0.275 s given when the intermediate 1 A singlet was first 
proposed. Redman et.al. deduced that a long lived 
electronic state accounted for a narrow 1.2 Hz resonance 
observed in a near degenerate 4- wave mixing experiment. 
They recognized that the narrow resonance could be as- 
sociated with a long lived spin state but their estimates 
of the spin-lattice relaxation times suggested otherwise. 
However, a spin-lattice relaxation time of 0.275 s is quite 
realistic for N-V centers in diamond and we consider that 



this was the correct interpretation. In this case their data 
can be explained without invoking a long lived singlet 
state and the information would be consistent with that 
presented here. 

The presence of a long lived singlet has also been 
adopted in accounting for the loss of emission from single 
centers as temperature is lowered 0, H| • With excitation 
it was considered that the singlet becomes populated but 
there is a thermally stimulated back transfer to the emit- 
ting level such that at room temperature there is little 
loss of emission. The back transfer decreases with low- 
ering temperature and the center can remain in the sin- 
glet level for a considerable time resulting in a drop of 
the average emission intensity. This decrease in intensity 
occurs for zero-phonon line excitation. However, if the 
explanation is correct there will be an equivalent loss of 
emission when the excitation is within the vibronic band. 
This is not the case. There is little change in N-V emis- 
sion intensity with a lowering of temperature (see Fig. 3). 
The more likely explanation for the loss of emission in 
the case of zero-phonon line excitation is spectral hole- 
burning. There are a range of processes (change of spin 
state, re-orientation of center or movement of charge in 
the neighborhood of the center) which can shift the ab- 
sorption frequency and cause a decrease of absorption 
for a laser held at fixed frequency. With a decrease in 
absorption there will be an equivalent decrease in emis- 
sion and the effect will become more pronounced as the 
temperature is lowered due to narrowing of the homoge- 
neous line width. The decrease in emission is, therefore, 
attributed to this process and not with populating the 
singlet level. Clearly no information about the energy of 
a singlet can be obtained from such experiments. l42| . 

It is often assumed that there is no decay from the sin- 
glet to the ground state 0, . Should this be the case 
the spin polarization would have to arises through the 
back transfer and such a process would be strongly tem- 
perature dependent. This is not what is observed. Spin 
polarization occurs from liquid helium temperatures to 
room temperatures and in the original measurements of 
spin polarization Loubser and van Wyk [T^ have shown 
that the polarization is maintained to 500 K. Another 
more important issue is to question how spin polariza- 
tion arises with the thermal back transfer process. No 
details have been presented as to how the spin polariza- 
tion occurs. This is in contrast with the present model 
where, rather than back transfer, there is decay from 
the singlet direct to the ground state. With the decay 
channels proposed one can readily account for the spin 
polarization. It can be concluded that the N-V center 
can be understood without back transfer playing a role. 

It is noted that the analysis of the photon statis- 
tics associated with emission from single N-V centers 
S [J S E3 has assumed that either there is signifi- 
cant singlet-triplet back transfer or a long lived singlet 
state (or both). Consequently the parameters reported 
are inconsistent with the values obtained here. Such mea- 
surements need to be re-analyzed using the present model 
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and consideration given to contributions associated with 
photo-ionization. 

Wrachtrup and co-workers at University of Stuttgart 
and at the National Academy of Sciences of Belarus have 
presented an exciting range of single site measurements 
includin g de monstra ting asp ects of quantum computing 
0, S H US IH IH IM El El. In the model adopted 
to interpret their results they have spin as a good quan- 
tum number and spin-orbit interaction is neglected. It 
is doubtful that spin-orbit can be totally neglected but 
it is small and there is a correspondence in the energy 
level schemes and selection rules between their work and 
that given here. There are also some similarities in the 
values of the parameters. For example, we have deter- 
mined that the rate constant for inter-system crossing 
from the excited x',y' states to the singlet s has a value 
of k x / s — k y / a — 0.39 x 1/t. This is in good correspon- 
dence with a value of 0.5 x 1/t given by Nizovtsev et.al. 
[34|. They have also proposed that there is much slower 
transfer from the 3 E z' state to the singlet and we agree 
with this conclusion. They have given a value of k 2 ' S = 
2.5 x 10 -4 x 1/r. In our model it is considered to be 
zero but a small value such as they have given does not 
change the behavior of the system. There is, therefore, 
reasonable agreement when considering the populating of 
the singlet level. The situation for transfer out of the sin- 
glet level is different and there is some confusion. There 
is no agreement when the singlet is taken to have the 
long 0.275 s lifetime and when there is no decay from 
this singlet to the ground state [34| ■ In this situation the 
spin polarization and transfer out of the singlet has to be 
through thermal back transfer plus optically driven pro- 
cesses. However, these processes are not consistent with 
the optically induced spin polarization being independent 
of temperature and linear |24l l43l | . As there are no such 
processes in our model there can be no comparison made 
with the parameters given. In other work [S HSl US E3 
direct transfer from the singlet to the ground state is 
indicated and mention made of a short singlet lifetime 
|45| . No parameters are given to enable a comparison but 
clearly there is a consistency with the model presented 
here. 

Jelezko et al 0] have observed a single sharp zero- 
phonon line in the excitation spectrum of single centers. 
This is a very significant observation as it is crucial for 
readout for many N-V quantum information processing 
applications. Detection requires there to be a transition 
that cycles without change of spin projection and their 
observation indicates that the z <-> z' transition can cycle 
^100,000 x before the change occurs |4^. This is in con- 
trast to that obtained by our model, where with present 
parameters, the cycling of the z <-> z' transition would be 
limited to ~50 x before a change of spin state. In the 
model the cycling is limited by the non-spin-conserving 
optical transitions and, as noted earlier, the strength of 
these transitions vary with strain. Strain varies the sep- 
aration of S z and (S x , S«) spin levels and consequently 
the degree of mixing via the non-axial spin-orbit interac- 



tion. This variation has been investigated experimentally 
by Santouri et al |2^. By studying small regions of an 
irradiated crystal they were able to obtain an inhomoge- 
neous line width (15 GHz) orders of magnitude narrower 
than obtained previously. Changing the spacial location 
gave spectra for different magnitudes of strain and the 
3 E splitting was resolved. Furthermore, clear hyperfine 
structure associated with the various optical transitions 
was obtained from two-laser hole burning measurements. 
Previous two-laser hole burning data indicated that the 
S z spin state was lowest in both the lower and upper 
component of the strain-split 3 E state. This is confirmed 
in the recent hole burning measurements and it is shown 
that this is the case with even small strain splittings of 
< 10 GHz. With such strain fields the order of the levels 
is, therefore, dominated by " spin-spin type" terms rather 
than by spin-orbit. The effect of this interaction at zero 
strain in Fig. 1(a) is to displace the S^, S y states with 
respect to the S z . This increases the spin separations in 
the upper branch and reduces the separation in the lower 
branch (or vice versa) and with increasing strain there 
will be a crossing of the spin levels in the lower energy 
branch (between Fig 1(a) and 1(b)). Thus, in the lower 
branch the levels are closer and give larger mixing of S z 
and S x , S y spin states. Transitions to the lower branch 
are, therefore, conducive to hole burning and electromag- 
netic induced transparency |29| as both require non-spin 
conserving transitions. Alternatively transitions to the 
upper level are more favorable for cyclic transitions and 
with low strain such transitions may account for the very 
cyclic transitions observed by Jelezko et. al. [9j. An 
understanding of these processes and their variability re- 
quire knowledge of the magnitude of the interactions as- 
sociated with the 3 E excited state. This has not been 
obtained in detail and remains an outstanding issue for a 
full understanding of the electronic structure of the N-V 
center. 



VIII. CONCLUSION 

In this work group theoretical considerations have been 
used to account for the electronic states of the N-V cen- 
ter, for the optical transitions and for the inter-system 
crossing. The account leads to a seven level model 
where the dynamics are dominated by four rate con- 
stants. These four rate constants are determined from 
independent experimental measurements and it is shown 
that with the values obtained the model gives plausible 
correspondence with additional optical measurements of 
ensembles. The comparison between theory and experi- 
ment is not fully quantitative but the agreement is suf- 
ficient to give confidence in the appropriateness of the 
model. A good physical understanding of the response 
of the center to optical excitation is obtained and the 
significance is that the model provides a basis for the 
development of strategies to target the remaining out- 
standing issues regarding the properties of the N-V cen- 
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tcr. The model also provides a sufficient understanding 
of the dynamics to allow for satisfactory development of 
many N-V applications. 
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